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1.  Introduction 

One  important  aspect  of  ocean  model  design  is  the 
choice  of  the  vertical  coordinate  system,  Traditional 
ocean  models  use  a  single  coordinate  type  to  represent  the 
vertical,  but  model  comparison  exercises  performed  in 
Europe  (Dynamics  of  North  Atlantic  MOdels  - 
DYNAMO)  (Willebrand  et  a!.,  2001)  and  in  the  United 
States  (Data  Assimilation  and  Model  Evaluation 
Experiment  -  DAMPE)  (Chassignet  et  al„  2000)  have 
shown  that  none  of  the  three  main  vertical  coordinates 
presently  in  use  (depth  [z-levels],  density  [isopycnal 
layers],  or  terrain-following  [cr-levels])  can  by  itself,  be 
optimal  everywhere  in  the  ocean.  The  HYbrid  Coordinate 
Ocean  Model  (HYCOM)  (Bleck,  2002)  is  configured  to 
combine  ail  three  of  these  vertical  coordinate  types.  It  is 
isopycnal  in  the  open,  stratified  ocean,  but  uses  the 
layered  continuity  equation  to  make  a  dynamically 
smooth  transition  to  a  terrain-following  coordinate  in 
shallow  coastal  regions,  and  to  z- level  coordinates  in  the 
mixed  layer  and/or  unstratified  seas.  The  hybrid 
coordinate  extends  the  geographic  range  of  applicability 
of  traditional  isopycnic  coordinate  circulation  models 
toward  shallow  coastal  seas  and  unstrati  fied  parts  of  the 
world  ocean.  It  maintains  the  significant  advantages  of  an 
isopycnal  model  in  stratified  regions  while  allowing  more 
vertical  resolution  near  the  surface  and  in  shallow  coastal 
areas,  hence  providing  a  better  representation  of  the  upper 
ocean  physics,  HYCOM  is  designed  to  provide  a  major 
advance  over  the  existing  operational  global  ocean 
Prediction  systems,  since  it  overcomes  design  limitations 
t*ie  present  systems  as  well  as  limitations  in  vertical 
^d  horizontal  resolution.  The  result  should  be  a  more 
streamlined  system  with  improved  performance  and  an 
Wended  range  of  applicability  (e.g„  the  present  systems 
re  seriously  limited  in  shallow  water  and  in  handling  the 
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transition  from  deep  to  shallow  water).  Figure  I 
illustrates  HYCOM ’s  ability  to  simulate  both  deep  and 
shajlow  currents  in  the  northern  South  China  Sea  (SCS). 

Global  HYCOM  with  1/12°  horizontal  resolution  at 
the  equator  (-7  km  at  mid-latitudes)  will  be  the  ocean 
model  component  of  an  eddy-resolving  operational 
nowcast/ forecast  system  scheduled  for  transition  to  the 
Naval  Oceanographic  Office  (NAVOCEANO)  at  the  end 
of  FY  07.  It  will  provide  nowcasts  and  forecasts  of  the 
three-dimensional  global  ocean  environment.  HYCOM 
will  be  coupled  to  the  Los  Alamos  CICE  sea-ice  model 
(Hunke  and  Lipscomb,  2004)  via  the  Earth  System 
Modeling  Framework  (ESMF)  (Hill  et  al.,  2004). 
Coupling  between  the  ocean  and  ice  models  will  more 
properly  account  for  the  momentum,  heat  and  salt  fluxes 
at  the  ocean/ice  interface.  The  final  component  of  the 
nowcast/forecast  system  is  the  NRL  Coupled  Ocean  Data 
Assimilation  (NCODA)  which  is  a  multivariate  optimal 
interpolation  scheme  that  will  be  used  to  assimilate 
surface  observations  from  satellites,  including  altimeter 
and  Multi-Channel  Sea  Surface  Temperature  (MCSST) 
data,  sea  ice  concentration  and  also  profile  data  such  as 
expendable  bathythermographs  (XBTs),  conductivity 
temperature  depth  (CTDs)  and  ARGO  floats  (Cummings, 
2006).  '  By  combining  these  observations  via  data 
assimilation  and  using  the  dynamical  interpolation  skill  of 
the  model,  the  three-dimensional  ocean  state  can  be 
accurately  nowcast  and  forecast.  Further  details  regarding 
all  data  assimilation  techniques  that  may  be  incorporated 
in  HYCOM  can  be  found  in  Chassignet  et  al.  (2006). 


271 


F'gtJr8  1:  *70p>  Mean  uPPer  ,ayer  currents  in  the  northern 
South  China  Sea  (SCS)  from  1/12'  global  HYCOM  using 
climatological  wind  and  thermal  forcing  from  the  ECMWF 
Note  the  across  shelf  flow  as  the  SCS  branch  off  the 
Kuroshio  feeds  the  northeastward-directed  SCS  Warm 
Current  along  the  Chinese  coast,  which  flows  counter  to  the 
prevailing  northeast  to  southwest  directed  monsoon  winds 
(Bottom)  Zonal  mean  velocity  along  117*E  in  the  upper  500  m 
with  blue  -  westward  and  yellow/orange  =  eastward.  HYCOM 
makes  a  smooth  transition  from  the  deep  to  shallow  water.  It 
®aPabl8  of  handling  the  sharp  topographic  rise  near 

Zl  N.  The  black  line  Indicates  the  location  of  the  vertical 

cross-section. 


2.  Model  Setup 


The  global  model  is  configured  on  a  Mercator  grid 
from  78°S  to  47°N,  while  north  of  this  latitude  an  Arctic 
dipole  patch  is  used  to  avoid  the  singularity  at  the  pole. 
The  current  1/12°  equatorial  resolution  translates  to  an 
array  size  of  4500*3298.  The  majority  of  the  experiments 
performed  this  year  used  32  hybrid  layers  in  the  vertical. 
Globa)  HYCOM  is  run  on  784  dedicated  processors  (98 
nodes)  on  the  NA VOCEANO  IBM  Power  4+  (kraken).  A 
typical  model  month  of  integration  uses  -22  wallclock 
hours  and  generates  525  Gb  of  uncompressed  output 
(-25 1  Gb  of  compressed  output), 

HYCOM  was  initialized  using  temperature  and 
salinity  from  the  1/4°  Generalized  Digital  Environmental 
Model  (GDEM3)  climatology.  So  far,  the  majority  of  the 
global  experiments  have  used  climatological  monthly 
mean  wind  and  thermal  forcing  constructed  from  the 
1.125°  European  Center  for  Medium-Range  Weather 
Forecasts  (ECMWF)  Re-Analysis  (ERA  1 5)  over  the 
1979—1993  time  frame.  In  addition,  6- hourly  variability 
Irom  the  ECMWF  operational  model  over  the  period 
September  1 994-September  1995  was  added  to  the 
climatological  wind  forcing  (but  not  the  thermal  forcing) 
to  add  the  higher  frequency  variability  needed  for  realistic 


simulation  of  the  surface  mixed  layer.  In  order  to  ke 
the  evaporation-precipitation  budget  on  track,  the  modH 
weakly  relaxes  to  the  Polar  Science  Center  Hydrographic 
Climatology  sea  surface  salinity  (Steele  et  al.,  200 1 
Global  HYCOM  does  not  currently  assimilate  anv 
oceanic  data,  but  that  capability  will  be  added  very  soon 
Assimilative  results  presented  below  come  from  a  Gulfnf 
Mexico  version  of  HYCOM, 

3.  Non  Assimilative  Global  HYCOM 
Evaluations 


Evaluation  and  validation  are  important  components 
of  an  ocean  model’s  developmental  process.  Comparing 
the  model  with  observational  data  may  highlight  both  its 
strengths  and  weaknesses.  Such  evaluations  may  guide 
the  developers  and  help  focus  (he  effort  to  where  it  is 
needed  most.  An  important  first  order  test  of  any  ocean 
model  is  proper  representation  of  the  large  scale 
circulation  and  eddy  structure.  This  is  illustrated  in 
Figure  2  that  compares  sea  surface  height  (SSH) 
variability  derived  from  satellite  altimeter  data  and  from 
1/12°  global  HYCOM.  The  model  is  able  to  capture 
much  of  the  expected  eddy  structure,  e.g.,  within  the 
western  boundary  current  regions  of  the  Gulf  Stream  and  | 
Kuroshio  Extension,  in  the  Gulf  of  Mexico,  along  the 
length  of  the  Antarctic  Circumpolar  Current  and  in  the 
Brazil-Malvinas  region.  Realistic  representation  of  the 
eddy  field  implies  the  model  is  accurately  simulating  the 
large-scale  flow  fields.  The  observed  SSH  variability  in 
the  tropica!  latitudes  is  not  reproduced  in  the  model 
because  it  was  am  with  climatological  wind  and  thermal 
forcing,  thus  excluding  the  interannual  variability 
associated  with  El  Nino. 


Figure  2.Sea  surface  height  (SSH)  variability  from  (top) 
observed  satellite  altimeter  data  based  on  Topex-Poseidon, 
ERS-1  and  ERS-2  over  the  period  October  1992  to  May  2005 
(derived  from  Collecte  Localisation  Satellite  [CLS],  France) 
and  (bottom)  simulated  from  1/12“  global  HYCOM  using 
ECMWF  climatological  wind  and  thermal  forcing 


! 
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Looking  below  the  surface,  global  HYCOM  is 
genera  lly  able  to  realistically  simulate  the  subsurface 
current  structure.  Comparisons  between  climatological 
observ  ed  and  simulated  currents  are  routinely  analyzed 
throughout  the  course  of  a  model  simulation.  Space 
limitations  prevent  a  complete  survey  of  these  results,  but 
HYCOM  realistically  simulates  the  surface  and 
subsurface  tropical  current  systems  (especially  in  the 
Pacific  Ocean),  the  Kuroshio  as  it  enters  Luzon  Strait  that 
connects  the  Pacific  Ocean  and  the  South  China  Sea 
(Figure  3),  the  downstream  flow  east  of  Taiwan  and  that 
entering  the  East  China  Sea,  and  the  Kuroshio  south  of 
Japan.  While  the  structure  of  the  simulated  upper  ocean 
velocity  field  in  the  Intra-Americas  Sea  is  consistent  with 
observations,  the  magnitudes  are  somewhat  weak  in 
global  HYCOM.  Efforts  are  currently  underway  to 
correct  this  deficiency. 


Figure  4.  Histograms  of  (left)  root  mean  square  difference 
and  (right)  correlation  coefficient  between  126  observed 
coastal  and  island  sea  level  stations  and  1/12*  global  HYCOM 
using  3-hourly  FNMOC  NOGAPS  wind  and  thermal  forcing  for 
2003.  A  30-day  running  average  is  applied  to  ail  time  series. 
The  median  rms  difference  is  4.4  cm  and  the  median 
correlation  coefficient  t$  .85. 


4,  Data  Assimilative  Progress  Using  NCODA 


Figure  3,  (Right)  Mean  upper  layer  currents  in  Luzon  Strait 
from  1/12*  global  HYCOM  using  climatological  wind  and 
thermal  forcing  from  the  ECMWF.  (Left  -  top)  Mean 
meridional  velocity  in  the  upper  300  m  based  on  a  ten-year 
composite  of  shipboard  acoustic  Doppler  current  profiler 
data  from  Liang  et  at.  (2003)  and  (left  -  bottom)  from  global 
HYCOM.  Both  observed  and  simulated  currents  indicate  a 
spilt  Kuroshio  with  a  stronger  western  than  eastern  branch 
and  weak  counterflow  in  the  lee  of  the  Babuyan  Island  chain. 
The  black  line  indicates  the  location  of  the  vertical  cross- 
section. 


One  simulation  was  performed  that  extended  a 
climato logically  forced  experiment  with  interannuai  wind 
and  thermal  forcing  from  the  2003-2005  3-hourly  Fleet 
Numerical  Meteorology  and  Oceanography  Center 
(FNMOC)  Navy  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS),  Interannually  varying 
forcing  allows  for  direct  contemporaneous  comparison 
between  simulated  and  observed  phenomena.  One  such 
validation  examines  coastal  and  island  sea  level 
variability.  Because  this  signal  is  a  deterministic  response 
of  the  wind-driven  circulation  at  many  locations,  a  non- 
assimilative  model  may  be  expected  to  perform 
reasonably  well.  This  is  corroborated  in  Figure  4  which 
shows  the  root  mean  square  (RMS)  difference  and 
correlation  coefficient  between  global  HYCOM  and  tide 
gauge  stations. 


Because  of  the  high  computational  requirements  for 
integration  of  1/12°  global  HYCOM,  development  and 
testing  of  the  NCODA  scheme  has  been  performed  using 
a  1/12°  Gulf  of  Mexico  version  of  HYCOM.  This  allows 
for  more  rapid  turnaround  and  was  performed  outside  our 
DoD  Challenge  project.  Scaling  it  up  to  the  global  model 
will  be  part  of  the  challenge  project. 

NCODA  is  an  oceanographic  version  of  the 
multivariate  optimum  interpolation  (MV 01)  technique 
widely  used  in  operational  atmospheric  forecasting 
systems.  The  ocean  analysis  variables  in  NCODA  are 
temperature,  salinity,  geopotential  (dynamic  height),  and 
velocity.  The  horizontal  correlations  are  multivariate  in 
geopotential  and  velocity,  thereby  permitting  adjustments 
to  the  mass  field  to  be  correlated  with  adjustments  to  the 
flow  field.  NCODA  assimilates  all  available  operational 
sources  of  ocean  observations.  Figure  5  shows  an 
example  of  the  Gulf  of  Mexico  HYCOM-NCODA 
cycling  system.  The  independent  overlaid  frontal 
locations  clearly  indicate  the  assimilation  is  keeping  the 
SSH  field  on  track.  Additional  experiments  assimilating 
simulated  “observational”  SSH  and  SST  data  starling 
from  an  incorrect  initial  state  demonstrates  that  the 
assimilation  reduces  the  subsurface  temperature,  salinity 
and  velocity  errors  after  subsequent  cycling  (not  shown). 
The  simulated  data  are  derived  from  a  control  run  used  as 
truth  in  verifying  convergence  of  the  data  assimilative  run 
toward  the  time  varying  SSH  and  3D  temperature,  salinity 
and  velocity  fields  of  the  control  run. 
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Figure  5,  Nowcast  SSH  at  four  dates  from  1/12*  Gulf  of 
Mexico  H'YCOM  using  the  NCODA  scheme.  The  HYCOM- 
ISCOOA  cycling  system  assimilates  SSH  at  altimeter  ground 
tracks,  MCSST  and  XBT  profile  data.  Overlaid  are  frontal 
locations  based  on  MCSST  data  and  produced  by 
NAVOCEANO.  Isolated  points  define  eddy  centers.  White 
(black)  lines  or  points  are  newer  (older)  than  four  days. 


5.  Plans 


By  the  end  of  F Y  06  we  expect  to  have  a  May  2001 
to  present  data-assimilative  hindcast  and  to  start  daily 
near  real-time  runs  that  mimic  the  expected  operational 
procedure.  In  FY  07,  we  plan  to  complete  a  1 993 -present 
"ocean  re-analysis”  by  running  a  hindcast  from  1993  up 
through  May  20GL  In  addition  we  will  complete  the 
1995-2007  NOGAPS-forced  non- assimilative  simulation 
and  perform  a  similar  simulation  from  1979-2006  with 
ECMWF  forcing.  These  interannual  simulations,  and  the 
data -assimilative  1993-present  hindcast,  are  an 
unprecedented  opportunity  to  study  regional 
oceanography  and  dynamics  at  high  resolution  in  the 
context  of  the  global  ocean. 
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